An Inverse Boundary Value Problem for the 
Magnetic Schrodinger Operator on a Half 

Space 

Valter Pohjola 
o ; 

D. September 6, 2012 

in 



University of Helsinki 
Department of Mathematics and Statistics 

Licentiate Thesis 
Advisors: Katya Krupchyk, Lassi Paivarinta 



Abstract: This licentiate thesis is concerned with an inverse boundary 
value problem for the magnetic Schrodinger equation in a half space, for 
potentials A G W^ p (W,R 3 ) and q e L? omp (W,C)- We prove that q and 
the curl of A are uniquely determined by the knowledge of the Dirichlet-to- 
Neumann map on parts of the boundary of the half space. The existence and 
uniqueness of the corresponding direct problem are also considered. 
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1 Introduction 



The main purpose of this thesis is to investigate an inverse problem for the 
magnetic Schrodinger operator in the half space geometry. The magnetic 
Schrodinger operator L^ q is defined by 

3 

L A , 9 :=^H9 J+ A 3 ) 2 + 9 (x). (1.1) 
j"=i 

We fix the half space by considering R 3 _ := {i 6 M 3 | 13 < 0}. We shall 
moreover assume that 



A G W^ p (Rl,R 3 ), and q G L~ mp (R 3 _,C), Imq < 0. (1.2) 



Here 



^com P ( R -' R3 ) := i A \w I A e W 1,0 °(K 3 ,M 3 ),supp(A) C M 3 compact} 
and similarly, we define 

L™ omp (W,C) ■= {q G L°°(R3,C) I supp(g) C W compact}. 

The direct problem, from which the inverse problem stems, is the Dirichlet 
problem 

(La „ — k 2 )u = inR 3 , 

I f (1.3) 

where k > is fixed and / G Hc/,m P (dM. 3 _). Furthermore, we will also require 
that the solution u should satisfy a boundary condition at infinity, which will 
be the Sommerfeld radiation condition 

lim \x\( -iku(x) ) = 0. (1.4) 



\x\->oo \ d\x 

Solutions satisfying this condition are called outgoing or radiating solutions. 
We will also occasionally use the term incoming solution. This refers to a 
solution of (ll.3p that satisfies (11.4)) . when the factor — ik is replaced by ik. 

The following result, which will be established in the first part of this 
work, gives the solvability of the direct problem (11.3)) . (I1.4p . 

Theorem 1.1. Let A G W^£ p (W_, M 3 ) and q G L™ mp (W_X) be such that 

Imq < 0. Then for any f G Hcllp(dRi), the Dirichlet problem (TO]) . (TO)) 
/ias a unique solution u G i?^ c 
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Here 

HIM) ■= U\w\ u ^ H L(^)}- 

Theorem II. II permits us to define the so called Dirichlet to Neumann map 
A A>q , (DN-map for short), A A>q : H^(dRl) ->■ H^(dRl) as 

/ >-> (d n + %A ■ n)u\ dR 3_, 

where u is the solution of the Dirichlet problem (11.31) . fll.4p and / is the 
boundary condition. Here n = (0, 0, 1) is the unit outer normal to the bound- 
ary dR 3 _. 

The inverse problem is then to determine if the DN-map uniquely deter- 
mines the potentials A and q in IR 3 _ . It turns out that the DN-map does not 
in general uniquely determine A. This is due to the gauge invariance of the 
DN-map, which was first noticed by [31]. It follows from the identities 

e~^L Aq e^ = L A+V ^ q , e~^k A ^ = A A+W ,„ (1-5) 

that A Aj(] = Aa+vv^ wnen ip G C 1 ' 1 (M 3 , 1R) compactly supported is such that 
^IdR 3 =0 (see Lemma ETTJ. This shows that A A>q cannot uniquely determine 
A. The DN-map does however carry enough information to determine VxA, 
which is the magnetic field in the context of electrodynamics. 

We shall use the notation Aj = (A/,i, Aj j2 , Aj t3 ) for the component func- 
tions, when considering a pair of magnetic potentials Aj, j = 1,2. We now 
state the main result of this work, which generalizes the corresponding results 
of [21], obtained in the case of the Schrodinger operator without a magnetic 
potential. 

Theorem 1.2. Let Aj e W^ p (W,R 3 ) and qj E L^ omp (W,C) be such that 
Imqj < 0, j = 1, 2. Denote by B an open ball in R 3 , containing the supports 
of Aj, and qj, j = 1,2. Let ri,r 2 C <9IR 3 be open sets such that 

(dR 3 . \B) DTj y£ 0, j = l,2. (1.6) 

Then if 

A Allfll (/)|r 1 = A^ iB (/)| riJ (1.7) 
for all f E Hclm P {dRl), supp(/) C T 2 , then 

V x Ai = V x A 2 and q x = q 2 in R 3 _. 

We would like to emphasize that in Theorem II .2[ the set Fx, where mea- 
surements are performed, and the set T 2 , where the data is supported, can 
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be taken arbitrarily small, provided that (11 .6p holds. The result of Theorem 
11.21 pertains therefore to inverse problems with partial data. Such problems 
are important from the point of view of applications, since in practice, per- 
forming measurements on the entire boundary could be impossible, due to 
limitations in resources or obstructions from obstacles. 

The first uniqueness result, in the context of inverse boundary value prob- 
lems for the magnetic Schrodinger operator on a bounded domain, was ob- 
tained by Sun in [31], under a smallness condition on A. Nakamura, Sun 
and Uhlmann proved the uniqueness without any smallness condition in [27] , 
assuming that A e C 2 . Tolmasky extended this result to C 1 magnetic poten- 
tials, in [35] and Panchenko to some less regular but small magnetic poten- 
tials in [22] . Salo proved uniqueness for Dini continuous magnetic potentials 
in [21]. The most recent result is given by Krupchyk and Uhlmann in [20] , 
where uniqueness is proved for L°° magnetic potentials. In all of these works, 
the inverse boundary value problem with full data was considered. 

In [10], Eskin and Ralston obtained a uniqueness result for the closely 
related inverse scattering problem, assuming the exponential decay of the 
potentials. The partial data problem in the magnetic case was considered by 
Dos Santos Ferreira, Kenig, Sjostrand and Uhlmann in [9] and by Chung in 

The inverse problem for the half space geometry, without a magnetic 
potential was examined by Cheney and Isaacson in [3]. The uniqueness 
for this problem in the case of compactly supported electric potentials was 
proved by Lassas, Cheney and Uhlmann in [21], assuming that the supports 
do not come close to the boundary of the half space. The result of Theorem 
11.21 is therefore already a generalization of the work [21], even in the absence 
of magnetic potentials. Li and Uhlmann proved uniqueness for the closely 
related infinite slab geometry with A = 0, in [2B]. Krupchyk, Lassas and 
Uhlmann did this for the magnetic case in [19]. In both of these works, the 
reflection argument of Isakov [16] played an important role. The uniqueness 
problem for the magnetic potentials in the slab and half space geometries 
has also been studied in a recent paper by Li [25]. The half space results in 
[25] differ from the ones given in this work, by concerning the more general 
matrix valued Schrodinger equation and by assuming C 6 regularity on the 
magnetic potential. 

The half space is perhaps the simplest example of an unbounded region 
with an unbounded boundary. It is of special interest in many applications, 
such as geophysics, ocean acoustics, and optical tomography, since it pro- 
vides a simple model for semi infinite geometries. We would like to mention 
that the magnetic Schrodinger equation is closely related to the diffusion ap- 
proximation of the photon transport equation, used in optical tomography 
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[3]. The half space geometry is also of interest in optical tomography, since 
in practice, the source-detector pairs are often located on the same interface 

The thesis is divided into two main parts. Section 2 gives a detailed ac- 
count of the solvability of the direct problem and provides a proof of Theorem 
11.11 Subsection 12. II develops some basic tools from scattering theory In Sub- 
section 12.21 we prove the existence and uniqueness for the direct scattering 
problem in all of 1R 3 , using the Lax-Phillips method. In Subsection 12.31 we 
extend this discussion to the half space case, using a reflection argument. 

Section 3 deals with the inverse problem and contains the proof of Theo- 
rem [L2j In Subsection 3.1 a central integral identity is derived. Subsection 
3.2 contains the construction of complex geometric optics solutions for mag- 
netic Schrodinger operators with Lipschitz continuous potentials. The proof 
of Theorem 11.21 is concluded in Subsections 3.3 and 3.4. 

1.1 Acknowledgements 

I wish to thank my advisor Katya Krupchyk for suggesting the topic of this 
thesis and for all the help while writing it. I also like to thank Mikko Salo 
for providing me with his notes on Carleman estimates for the Magnetic 
Schrodinger operator. 
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2 The direct problem 



The purpose of this section is to investigate the well-posedness of the bound- 
ary value problem ( II. 3p . (II .4]) for the magnetic Schrodinger operator in the 
half space, and to establish Theorem 11.11 

We prove existence and uniqueness results for this problem, but we will 
not concern ourselves with showing that solutions depend continuously on 
boundary data. This is enough to guarantee that the questions of the unique- 
ness of the inverse problem are sensible. 



2.1 The free space outgoing Green function 

The aim of this subsection is to introduce the outgoing Green function for 
the Helmholtz equation and develop some basic notions used in scattering 
theory that we need to attack the direct problem. For a more in depth 
exposition, see [7]. Having constructed the outgoing Green function, we use 
it to investigate the corresponding resolvent operator and the asymptotics of 
solutions. 

A function G is generally speaking a Green function for the Helmholtz 
equation if it solves the problem 

(-A x -k 2 )G(x,y) = 6(x-y), yeQ, (2.1) 

in some region Q and satisfies some specific boundary condition on dQ. We 
are interested in the case Q = M 3 with the boundary condition being the 
Sommerfeld radiation condition (II. 4p . 

In the next proposition we construct a specific Green function called the 
outgoing free space Green function denoted by G , which satisfies (12. ip and 
the Sommerfeld condition. 

Proposition 2.1. Let k > 0. The function 

ik\x-y\ 

G (x,y) = — (2.2) 

4ir\x — y\ 

is a free space outgoing Green function. 

Proof. We start by assuming y = in (12. ip and thinking of Go a s a function 
of x only. Instead of considering the operator —A — k 2 directly, we first look 
for outgoing Green's functions G e for the operators —A — k 2 —ie, where e > 0. 
By inserting ie into (12.11) and taking the Fourier transform we get that 

de - 1 

£ 2 — k 2 — ie 
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Since this is a locally integrable function, decaying at infinity, we may take 
the inverse Fourier transform. Let us therefore proceed by calculating the 
inverse Fourier transform, 



i p2iv pit poo &i£,(rfl>6)-x 

i. / / / o 9 

-r sin 



(2vr) 3 7 J J r 2 -k 2 - le 

The Fourier transform of a spherically symmetric function is spherically sym- 
metric. We need therefore to calculate this integral only for e.g. x = (0, 0, R), 
R = \x\. By doing this and abbreviating a = k 2 + ie we get 



o Jo 



co girR cos 9 j 

r 2 sin QdQdr 



r 2 — a iR J 

J roo ^—irR, gir_R 



' 7r ^irRcosO 

— rdr 

r 2 — a 



iR J r 2 — a 

1 f°° e irR 



rdr 



iR .Loo r 2 - a 



rdr. 



We can evaluate this last integral with the method of residues, if we tem- 
porarily replace r with a complex variable z. We choose as a contour, an 
origin centered half circle C s that is in the upper half plane and goes along 
the real axis. The parameter s stands for radius of the circle. By choosing 
the branch of the complex square root in the upper half plane and a large 
enough s, we achieve that the pole y/a will lie inside the contour. A residue 
integration gives us then 

rdr = lim <I> i=rzdz 

oo 



r 2 - a s->co J Cs (z - y/a) (z + y/a) 

2iri Res ( -=-. -=-, y/a 

z - y/a){z + y/a) 



By going back to (12.31) we see that 



1 \ JL\/k 2 +ieR 

^H- 2 — r2 — - )(x)= — — , R=\x\ 

X 2 - k 2 -iej AtxR 



In other words, G e (x) = e l ^ k ' 2+le \ x \ / (47r|x|). 
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We will need to use the theory of distributions to extend the above argu- 
ment to the case e = 0. We define the tempered distribution (£ 2 — k 2 — iO)^ 1 , 
by taking the limit 

((?-k?-iO)-\ <p )= lim {(e-k 2 -ie)-\ V ) (2.4) 

for tp G 5(IR 3 ), where the distribution pairing on the right is given by an inte- 
gral. One can check that this definition makes sense, using similar arguments 
as for the principal value distribution. We define 

Go -^^{(e-k'-ior 1 ). 

The first part of the proof shows that 

(F- 1 {(e - k 2 - zO)" 1 ) , ip) = lim ( T- x G t , <p) 

/ e iy/k 2 +ie\x\ 
; — —. (p(x)dx 
Atx\x\ 



e lk \x\ 
47r|x| 



<p(x)dx, 



so that G (x) = e ik ^/(A7r\x\). 

To check condition (12.11) . we first note that G e is a Green function, and 
hence 

lim ((-A-k 2 -ie)G t ,p) = (6,<p). (2.5) 
On the other hand we have 

lim ((-A - k 2 - ie)G e , if) = lim (((-A - k 2 )G e , (p) + ( - ieG ei ip)) 

<e->0+ v ' e->0+ V v ' v 'J 

= (G , (-A - k 2 )ip) + lim < - ieG e , ip). 

For the last term we have —ie{G e ,p) — > 0, as e — > 0. Combining the above 
with (12.51) gives 

((-A-k 2 )G ,p) = (5,p), 

so that (|2IIJ holds. 

A direct calculation shows that Go satisfies the Sommerfeld radiation 
condition (11.41) . Finally let us set Gq(x, y) = Gq(x — y). □ 
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Remark 2.2. Notice that G £ L 2 oc (IR 3 ). Smce G satisfies (12. ip m £/ie 
sense 0/ distributions, we have for ip £ L^ omp (IR 3 ) 

(-A - k 2 )(G * (p) = 5 * ip = if. 

Remark 12.21 motivates the notation, 

(-A - k 2 - iO)~VO) : = G * <p, 

where — iO marks the fact that we took the positive sided limit in (12.41) . 
This is important since the positive limit guarantees that Go and Go* if are 
outgoing, i.e. satisfy the Sommerfeld condition. We would have ended up 
with an incoming solution, had we taken the negative sided limit. 

The above operator is in fact a limiting value of the resolvent of the 
Laplace operator. This operator is not continuous on L 2 (IR 3 ), when Im k — 0. 
It is however continuous between certain weighted spaces (see e.g. [2]). 

The next lemma gives a more modest continuity result, which will be 
used later. 

Lemma 2.3. The operator 

(-A - k 2 - : L 2 comp (R 3 ) L 2 oc (R 3 ) 

is continuous. 

Proof. Let us write T := (—A — k 2 — iO)^ 1 . Let B r be a ball centered at 
the origin of radius r > and \B r stand for the characteristic function of 
the ball B T . It suffices to show that the operator Xb t Txb t is continuous on 
L 2 (R 3 ). We have, for x £ B r , 

TxbMx) = / G {x-y)tp{y)dy, up £ L 2 (R 3 ). 

Since Go £ L 2 (B r x B r ) we get by applying the Cauchy-Schwarz inequality, 
that 

/ \TXB r v(x)\ 2 dx < I ( / G (x-y)<p(y)dy) dx 

J B r J B r V J B r J 

< / \G (x - y)\ 2 dy \\ip\\ 2 2 dx 

«/ Bj- «/ Bj- 

<c\M 2 2 . 

Hence \\xB r TxB r ^h < C\\tp\\ 2 . 

□ 
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Next we use the closed graph theorem and elliptic regularity to extend 
the above result to case where the range is the Sobolev space Hf 1 



loc\ 



Corollary 2.4. The operator (-A - k 2 - iO)- 1 : L 2 comp {R 3 ) H 2 oc (R 3 ) %s 
continuous. 

Proof. Let us write T := (—A — k 2 — z'O) -1 . By elliptic regularity, we have 

TL 2 comp (R 3 ) C Hf oc (R 3 )- 

The claim follows from the closed graph theorem, once we show that the 
operator T : L 2 (B) -»■ H 2 (B) is closed. Here B is an open ball in R 3 . To 
that end, let f n G L 2 (B) be such that f n — > f in L 2 (B) and Tf n — >■ ^ in 
F 2 (S). By Lemma 1231 we know that Tf n ->■ T/ in L 2 (B). Hence, 5 = T/. 
This completes the proof. □ 

We now begin investigating the asymptotics of radiating solutions to the 
Helmholtz equation. First we look at the asymptotics of Go- 

Lemma 2.5. The outgoing Green function has the following asymptotics, 

ik(\x\-(x-y)/\x\) / 1 \ 

G ^-"» = 4,|x| +0 {w} (2 - 6) 

as \x\ — > 00, uniform for y in a bounded set. It can be differentiated with 
respect to y. 

Proof. We will derive the asymptotic expression estimating the nominator 
and the reciprocal of the denominator of Go separately. 
We start by writing 



\x - y\ 




By looking at the Taylor series we see that \/l + x = 1 + x/2 + 0(x 2 ), as 
x —7- 0. Applying this to the above gives 

1 1 11 x ~y , n( l \ 

\x — y\ = \x\ 1— r- + CM - — : , 

\x\ \\ X \J 

for |x| — > 00. It follows that 

1 ' - "I W (2.7) 



\x — y\ \x\ \\x 
12 



as \x\ — > oo. 

For the nominator of Gq, let f(x) = 0(1/ \x\) be the function, for which 
ik\x — y\ — ik\x\ — ikx ■ y/\x\ + f(x). Then 

exp(ik\x — y\) = exp NA;|x| — ik -^ exp(/(x)) 

= exp I ik\x\ — ik 1(1 + i\ x ) 



x \ V J v ' 2! 



which gives the expression 



x m y 

exp(ik\x — y\) = exp ( ik\x\ — ik- — - ) + O 



x\ J \\x 

as \x\ — > oo. Multiplying this with the asymptotic expression ( 12. 7p gives the 
asymptotic expression for Gq. □ 

The next small lemma shows that the L 2 norm over a sphere of an out- 
going solution to the Helmholtz equation stays bounded as the radius of the 
sphere grows. Note that this applies also to incoming solutions. This will be 
used later. The lemma is mainly needed, since we use a weak form of the 
Sommerfeld radiation condition (jl.4p . 

Lemma 2.6. Let u G Hf oc {E?) be an outgoing or incoming solution to the 
Helmholtz equation (—A — k 2 )u = in M 3 \ B r , where B r = {x | |x| < r}. 
Then 

\\u\\i?(8B.) = 0(l), \\d n u\\ L 2 (dBs) = 0(1), s -> oo. (2.8) 

Proof. Pick a ball B s D 5 r . Assume that tt is outgoing, i.e. it satisfies ( II. 4p . 
Multiplying this condition with its complex conjugate and integrating over 
B s gives 



/ (k 2 \u\ 2 + \d n u\ 2 + 2klm(ud n u))dS = \d n u - iku\ 2 dS -> 0, 

J\x\=s J \x\=s 

(2.9) 



as s — > oo. Here n is the unit outer normal to d B s . Using Green's formulas 
we have that 

/ (ud n u — ud n u)dS = I (uAu — uAu)dx = 0, 

JdB s UdBr JB s \B r 
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and therefore, 



(ud n u — ud n u)dS = — (ud n u — ud n u)dS 

dB r JdB s 



= —2i j lm(ud n u)dS. 

JdB s 

It follows that the last expression does not depend on s. This shows that the 
two first terms in (12. 9p . which are both positive, are bounded, and particu- 
larly we see that (12.81) holds. 

The incoming case can be reduced to the outgoing case as follows. Firstly 
if u is an incoming solution of the Helmholtz equation, then u is an outgoing 
solution. Moreover ||m||l2 = ||m||l 2 - The incoming case follows therefore from 
the outgoing case. □ 

The following lemma gives a boundary integral representation for radiat- 
ing solutions of the Helmholtz equation in exterior regions. 

Lemma 2.7. Let u G Hf oc {E?) be an outgoing solution to the Helmholtz 
equation (—A — k 2 )u = in M 3 \ B r , where B r = {x E M 3 \ \x\ < r}. Then 



u(x) = / (d ny G (x,y)u(y) - G (x,y)d ny u(y))dS(y), 

JdB ri 

for x6K 3 \ B ri and r±> r. 

Proof. Let r < r\ < r-i and Xq G M 3 be an arbitrary point in Q := B T2 \ B ri . 
Applying the Green formula to Q e := B r2 \ (B ri U B t (x )), e > sufficiently 
small, we get 

= / (— A — k 2 )Go(x 1 xq)u{x) — Gq{x, xq)(— A — k 2 )u(x)dx 

(2.10) 

(- d n G (x, x q )u(x) + G (x, x ) d n u(x))dS(x). 



We have 

Gq(x, xq) d n u(x)dS(x] 



ldB c {x ) 

Consider next 



< / < 0(e) 

8B t {x ) 47re 



d n Go(x, xo)u(x)dS(x) = — / — u(x)dS(x) 

dB c (x ) JdB e {x ) 4?re 



+ / - — ^u(x)dS(x). 
JdB £ (x a ) 47re 2 
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It follows that 



I 



ike 



ike 



dB e (x ) 



Aire 



■u(x)dS(x) 



0(e) 



as e — > 0. Using the fact that u[x) = u(x ) + 0(e), we get 



-u{x)dS{x) = e ike u(x ) + O(e). 



dB ( (x ) 



4ire 



Letting e — > in (12.10)) . we obtain that 



u(x ) 



d n Gq(x, xq)u(x) + Gq(x, x ) d n u(x))dS(x). 



a n 



The next step is to show that 



I :-- 



{— d n G (x, xq)u[x) + G (x, x ) d n u(x))dS(x) 0, 



as T2 — > oo. By adding and subtracting ikGou, I can be written as 



+ 



dB r 



dB r 



(d n Gq(x, xq) — ikGo(x, xo))u(x)dS(x) 



Go(x,xo)(d n ii(x) — iku(x))dS(x) 



We show that the first term goes to zero as r 2 — > oo. The second term can be 
estimated in the same way. Because of (12. 8p we can use the Cauchy-Schwarz 
inequality and write 



/ 

J dB r 



u(d n G — ikG )dS 



< 



[ \u\ 2 dS [ \d n G - ikG \ 2 dS 0, 



as r 2 — > oo. Here we have used that \d n Go — ikGo\ 2 = o(r 2 2 ), valid because 
of Qj). 

□ 

The preceding lemmas allow us to prove the main result on the asymp- 
totics of outgoing solutions to the Helmholtz equation. 

Proposition 2.8. Let u G Hf Q JJsL z ) be an outgoing solution to the Helmholtz 
equation (—A — k 2 )u = in 1R 3 \ B r , where 5 r = {x G I 3 |aj| < r}. Then 
there exists a G L 2 (S 2 ) such that 

p ik\x\ 



U[X) 



\x\ 



-a(x) + O 



1 



x\ 



x := x/\x\ G S 



(2.11) 



uniformly in all directions as x —> oo. 
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Proof. An application of Lemma 12.51 allows us to write 



£ ik\x\ 



G (x,y) = ——a^x.y) + fi(x,y), 



where f\ = 0(l/|x| 2 ) uniformly in all directions as \x\ — > oo. Next, a straight- 
forward computation using Lemma 12.51 shows that 



e ik\x\ 



d n G (x, y) = — r a 2 (x, y) + f 2 (x, y), 
\x\ 

where f 2 = 0(l/|x| 2 ) uniformly in all directions as |a;| — > oo. By the repre- 
sentation of Lemma 12.71 we have 

r / e ik\x\ \ 

u ( x ) = / -r-r a 2(^, V) + /2O, y) u(y)dS(y) 



dB ri 



\x\ 



dB ri 



e ik\x\ 



ai(x,y) + fi(x,y) )d n u(y)dS(y) 



with T\ > r. Let us split this into four separate integrations corresponding to 
the individual terms. The terms involving a\ and a 2 are clearly of the form 
of the first term on the right side of the equation (12.111) . The two remaining 
terms give the contribution, 



{h{x,y)u{y) - fi(x,y)d n u(y))dS(y). 



dB r 



Since fi, f 2 = 0{l/\x\ 2 ) uniformly in all directions as |x| — > 00, we conclude 
that the integral above is 0(l/\x\ 2 ), which proves the claim. □ 

The function a in (12.111) is called the far field pattern or scattering am- 
plitude and is of central interest in scattering theory. 



2.2 The magnetic Schrodinger equation in M 3 

In this subsection we begin proving existence and uniqueness for problem 
(11.31) . (II ,4p by first considering the partial differential equation in the whole of 
M 3 . The main tool will be the Lax-Phillips method from scattering theory, see 
[22] and [15]. We will assume that the potentials and sources are compactly 
supported in IR 3 . More precisely, we assume that 

A G W^ p (R 3 ,R 3 ), q G L^ mp (M 3 ,C), Imq < 0, and / G L c 2 omp (M 3 ). 

(2.12) 
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Our aim is to find an outgoing solution u to 

(L Ag -k 2 )u =/, inM 3 (2.13) 

and to show that it is unique. 

We begin by recalling what is meant by a weak solution. We call u G 
#/ oc (R 3 ) a weak solution to (IZTBl if for every v G C£°(M 3 ) 

/ (Vtt • Vv - 2iA ■ (Vit)iJ + P k ^ A uv)dx = / fvdx, (2.14) 

where Pk, q ,A '■= — ?V ■ A + A 2 + q — k 2 . We state explicitly the following 
regularity result for weak solutions. 

Lemma 2.9. Lei u be a weak solution to (I2.13p . Then u G Hf oc {M?). 
Proof. Since it G if/ oc (lR 3 ) and / G L 2 (IR 3 ), we have 

-Am = 2% A ■ S7u + P KqA u + f G ^L(^ 3 )- 
By elliptic regularity, see [12], we conclude that u G if ; 2 oc (lR 3 ). □ 

The asymptotics given by Proposition [231 combined with Rellich's lemma 
and the unique continuation principle, see the appendix, give the uniqueness 
for outgoing solutions to fl 2 . 1 3 [) . This is the content of the following theorem. 

Theorem 2.10. Assume that A andq satisfy f)2.12p . Ifu G if/ oc (M 3 ) satisfies 
the Sommerfeld condition (II. 4p and solves 

(L Ai g - fc 2 )it = in M 3 , 

where k > 0, then u = 0. 

Proof. Let B r := 5(0, r). Denote the L 2 (B r ) inner product by (•,•)• ft 
follows from the Green's formula of Lemma 14. 1| in the appendix that 

{{L A>i - k 2 )u,u) - (u, (L Aj0 - k 2 )u) - (u,qu) 
= (it, (d n +iA ■ n)u) L 2 (dBr) - ((d n +iA ■ n)u, u) L 2 (dBr) . 

The first term on the left side vanishes and the second term is (it, qu). Notice 
also that the vector field A vanishes along d B r . The above equation reduces 
thus to 

to/ «=to/ 9M ^<0. 
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Using the asymptotic expansions in Proposition 12.81 we get 

p -ik\x\ / 1 \ \ / ih P ik \ x \ f 1 

Im / [a e -^ + (A-A\ (a-^-^ + O A) )dS 



'dB r V l x l \l x l / / V \ x \ \\ x \ 

= Im J (^a{e,ip)\ 2 ^ + o(^^jr 2 sinO d6 dip. 
By taking the limit as r — > oo, we obtain that 

£;|a(#,^)| 2 sm«V < 0, 

l*l=i 

and hence the far-field pattern a vanishes identically. 

An application of Rellich's lemma, see Proposition 14.81 in the appendix, 
allows us to conclude that u = outside B r . The unique continuation prin- 
ciple, see Theorem 14.61 in the appendix, implies that u = in M 3 . □ 

We now proceed to proving the existence of outgoing solutions to (I2.13p . 
In doing so, we shall first establish that the Dirichlet realization of L^ q on a 
ball has a discrete spectrum. Showing this is complicated by the presence of 
the imaginary part of q, which makes LA, q non-self-adjoint. We will use the 
notion of relative compactness for operators to resolve this issue. 

Definition 2.11. Let B be a Banach space and let T be a closed densely 
defined linear operator on B such that Spec(T) ^ C. Assume that A is a 
linear operator on B such that D(T) C D(A). We say that A is relatively 
compact with respect to T if for any sequence {u n } C D{T), such that both 
{u n } and {Tu n } are bounded, {Au n } has a convergent subsequence. 

Lemma 2.12. Let B be a Banach space, let T be a closed densely defined 
linear operator on B, and assume that X Spec(T). Then A is relatively 
compact with respect to T if and only if A{\ — T) _1 is compact. 

Proof. Suppose that we have a sequence {u n } C D(T) such that {u n } and 
{Tu n } are bounded. Then there is a constant M such that 

|| (A — T)w n || < |A|||w n || + ||Tit n || < M < oo, 

for all n, so that {(A — T)u n } is bounded. Since A(X — T) _1 is compact, the 
sequence 

{A(\-T)-\\-T)u n } = {Au n }, 

has a convergent subsequence. 

A similar deduction shows the opposite direction of the lemma. □ 
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We will now show that the Dirichlet realization of L A , q on a ball has a 
discrete spectrum. 

Lemma 2.13. Assume that A and q satisfy (I2.12p and that 

supp(g), supp(A) C B r = {x E R 3 | \x\ < r}. 

The operator L Aq : L 2 (B T ) — > L 2 (B r ), equipped with the domain H 2 (B r ) n 
Hl(B r ), is closed and has discrete spectrum. 

Proof. We let L be the operator —A on L 2 (B r ), equipped with the domain 
H 2 (B r ) n Hq(B t .) =: D(Lq). We know that L is self-adjoint with discrete 
spectrum, see [T2"] . 

We will show that LA, q is a relatively compact perturbation of Lq. Ac- 
cording to Theorem 11.2.6 in [8], the essential spectrum is preserved under 
relatively compact perturbations, as is closedness. Thus L Aq will have no 
essential spectrum, since L has none. By writing 

L Aiq = -A - 2iA ■ V + p, p= -iV • A + A 2 + q e L°°(B r ), 

we see that our task is reduced to showing that —HA ■ V + p is relatively 
compact with respect to Lq. 

Assume that A ^ Spec(L ). By the criterion of Lemma [2.12} the operator 
—HA ■ V + p is relatively compact with respect to L if and only if 

(-2iA- V + p)(\I-L )- 1 

is compact on L 2 (B r ). We split this as 

-2iA ■ V(A7 - Lq)" 1 + p(XI - L y\ 

and show that both of the resulting operators are compact on L 2 (B r ). The 
resolvent operator (XI — Lo)" 1 is continuous: L 2 (B r ) H 2 (B r ) n H%(B r ). 
The latter space is however compactly imbedded into L 2 (B r ), by the Rellich- 
Kondrachov theorem. If we view p as a multiplication operator, then it 
is continuous: L 2 (B r ) — > L 2 (B r ). This shows that the second operator is 
compact. 

The operator 2iA ■ V is continuous from H 2 (B r ) — > H l (B r ), since A is 
Lipschitz. The latter space is however compactly imbedded into L 2 (B r ). The 
first operator is therefore also compact. 

□ 

We are now ready to prove the existence of outgoing solutions to (I2.13p . 
using the Lax-Phillips method. 
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Theorem 2.14. Assume that A and q satisfy (|2.12p . Let k > 0. Then 
for any f G L^ (M. 3 ), there exists u G Hf oc (MA), satisfying the Sommerfeld 
radiation condition, that solves 

(L A , g - k 2 )u = f in M 3 . 

Proof. Let B r = {x G IR 3 | |x| < r} be such that supp(g), supp(A) C B r . Let 
s > r be such that supp(/) C B s . We pick a function ip G C^°(M 3 , [0, 1]) 
such that <p = 1 on B r and supp((/3) C B s . 

Let A G C, Im A 7^ 0, be such that A avoids the spectrum of the Dirichlet 
realization of La q on the ball B s . The existence of such A is guaranteed by 
Lemma 12.131 

We begin by looking for a solution u of the form 

u = pw + (1 — <p)v. (2-15) 
Here w G H 2 (B S ) fl Hq(B s ) is the unique solution to the Dirichlet problem 

{L A , q - X)w = g in B s , 

w\ 9Bs = 0, 

where g G L 2 (IR 3 ), with supp(g) C B s . And v is the unique outgoing solution 
of the equation, 

(-A-k 2 )v = g inM 3 . (2.16) 



Remark 12.21 gives an explicit formula for v, and according to Corollary 12.31 
we know that v G i^ 2 oc (IR 3 ). 

Inserting u into the original equation will result in an operator equation 
for the unknown function g. Abbreviating L := {L A ,q — k 2 ), we have^ 

Lu = L(tpw) + L((l — ip)v) 

= [L, <p]w + ipLw + [L, (1 — ip)}v + (1 — p)Lv 

= [L, ip]w + (f(L A , q - X)w + cp(X - k 2 )w + [L, (1 - ip)]v + (1 - (p)g 

= [L, cp]w + (pg + (f(X - k 2 )w + [L, (1 - <p)]v + (1 - (f)g 

Noting that the commutator of k 2 and <p is zero, we get from the above, 

Lu = [L A>q , (p]w + (/?(A — k 2 )w + [L A>q , (1 — (p)]v + g 
= [L A , q , <p]w + <p(\ - k 2 )w - [L Aiq , <p]v + g. 



1 Here we use the following bracket notation for the commutator operator [A, B] := 
AB - BA. 
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By setting Tg := <p(X — k 2 )w + [L A ,q, — v ), we see that g is to satisfy 
the operator equation 

(I + T)g = f. (2.17) 

Our problem of finding a solution of the special form (I2.15P is thus reduced 
to showing that the equation (12.1 7p has a solution. 

Our aim is to use the Fredholm theory. In order to do this we need to show 
that T : L 2 (B S ) — > L 2 (B S ) is compact. Notice that the first term of T is y?(A — 
k 2 )(L Atq - \y 1 g. But since the resolvent (La* ~ A)" 1 : L 2 (B S ) -> H 2 (B S ) is 
bounded and since the inclusion map: H 2 (B S ) — > L 2 (B S ) is compact, we see 
that the first term of T is compact on L 2 (B S ). 

For the second term of T, we first note that (—A — k 2 — iO)" 1 : L 2 (B S ) — > 
H 2 (B S ) is bounded by Corollary 12.41 Next note that supp([LA,g, tp]v) C B s , 
since supp(y9) C B s . The first order operator [L A , q , is explicitly given by 

[L A , q , <p] = -A(^ - 2V<^ ■ V - 2iA ■ Vcp. 

It also has W 1,OQ coefficients, so that it maps H 2 (B S ) to H l (B s ) continuously. 
Now i/ 1 (i? s ) is compactly embedded in L 2 (B S ) and hence we see that T is 
compact on L 2 (B S ). 

According to the Fredholm theory, we need only to show that / + T is 
injective, to have the surjectivity and thus a solution to (I2.17p . for a given /. 

Assume that (J + T)g = 0. Theorem 12.101 gives that 

u = ipw + (1 -<p)v = inM 3 . (2.18) 

We want to show that g = 0. This follows if w = 0, which we will prove 
next. First note that (12.181) gives (fw = (</? — l)v , so that 

w = when tp = 1. (2.19) 

In particular, w = on 5 r . This gives that 

# = (L A>g - A)w = (-A - \)w 

in B s , because supp(A), supp(g) C B r . We also have that 

(~A-k 2 )v = g, 

in B s . Subtracting the last two equations and using that v = <p(v — w), we 
get 

(-A -£;> + (- A- \)(-w) = (-A-X)(v-w)-k 2 v + Xv 

= (-A-X)(v-w)-(k 2 -X)(p(v-w) 
= 0. 
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Set r :—v — w. Now (—A — A)r — {k 2 — \)<pr = in B s . Multiplying this by 
r, using integration by parts and noting that r\g Bs = — w\qb 3 = 0, because 
of fl2~T8j) we get 

/ |Vr| 2 - \\r\ 2 dx = I (k 2 - \)^\r\ 2 dx. 

J B s J B s 

Taking the imaginary part gives 

ImA / (1 - (p)\r\ 2 dx = 0. 
Jb s 

From this we see that r = in the region where tp ^ 1. It follows that 
v = tpr = when tp ^ 1. Hence, 

w = v — r = when tp ^ 1. (2.20) 

Combining fl2TT9|) and f l2T20|) . we see that w = in B s . 

□ 

In summary, we see that Theorem 12.101 and Theorem 12.141 show the exis- 
tence and uniqueness of outgoing solutions for the problem f)2.13p . 



2.3 The magnetic Schrodinger equation in a half space 

The main objective of this subsection is to extend the existence and unique- 
ness results of the previous subsection to the half space case. 
Let k > and let 

AeW^ p (W,R 3 ), geLl p (#,C), lmg<0. (2.21) 

Given / e L 2 comp (W_X) ■= {/ G L 2 (W,C) | supp(/) C Incompact}, we 
first prove the existence and uniqueness of an outgoing solution to the fol- 
lowing problem, 

(L Aa -k 2 )u= f in R 3 , 
I n ( 2 - 22 ) 

We will reduce this problem to the case of M 3 , by using an extension argument 
and then use Theorems 12.101 and 12.141 The extension argument relies on the 
following Lemma, which will also be of importance later (see also Theorem 
1.3.3 in P3|). 
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Lemma 2.15. Let v G C 0,1 (IR 2 ), with compact support. Then there is a 
i/) € C 1 ' 1 (IR 3 ), with compact support for which 

t/j{x,0)=0 and d 3 ip(x, 0) = v(x), (2.23) 

for x e M 2 . 

Proof. Let ip be the usual mollifier function in M 2 , with (p e C^°(M 2 ), </? > 
and J cp — 1. Set = l/t 2 (p(x/t). We define u(x,t) := (v * ip t )(x), for 

t ^ and m(x, 0) := t>(x) = (t> * 5)(x). More explicitly 

M (M) = i / v{y)<p(^-—^)dy= v(x - ty)tp(y)dy, (2.24) 

1 Jl2 V t / ,/ R 2 

for t ^ 0. From the right hand side we see that u is Lipschitz in (x,t), 
because v is Lipschitz. We define ^ as 

V(x,t):=tu(x,t). (2.25) 

Notice that ^> satisfies the first condition in f)2.23p . 

Next we show that the partial derivatives of \1/ are Lipschitz. When t ^ 
we have, using f)2.24p that 

d Xi (tu) = t J v{y){di if) ^dy 

v(x - ty) di(p(y)dy. 

It follows that this identity also holds when t = 0. The right hand side of this 
identity is easily seen to be Lipschitz in both x and t, since v is Lipschitz. It 
follows that d Xi (tu) is Lipschitz in M 3 . 

The next step is to show that dt \& is Lipschitz. To see that dt $ = dt(tu) 
is continuous, we compute the derivative at t — 

_ , , hu(x, h) — u(x, 0) . 
d t (tu)\ t=0 = lim 1 ; = u(x,0) = 

Notice that this shows that \l/ also satisfies the second condition in f)2.23p . 
And further that d t (tu) is Lipschitz in x when t = 0. 

For t 7^ 0, observe firstly that d t (tu) = u + td t u. Thus we need only to 
check that the later term is Lipschitz. We write using (12.241) 

d t u = - [ v{x-ty)(Vip(y)-^ + <p{y)l)dy. (2.26) 
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This gives that 

td t u = - / v(x -ty)(V<f(y) ■ y + 2<p(y))dy, 

which is easily seen to be Lipschitz in both x and t, because v is Lipschitz. 

To obtain ip we pick a % G Cjf (M 3 ), s.t. x\ supp(t') = 1. Then ip := \^ G 
C 1 ' 1 ^ 3 ), X*\t=o = and d t ( X *)\t=o = (x <9 t *)|t=o = v. 

□ 

We use the above Lemma to show that, it is sufficient to consider existence 
and uniqueness in problem (12.221) . for potentials A for which d 3 A\ X3= o = 
0. To see this notice that Lemma 12.151 guarantees the existence of a ip G 
C 1,1 (M 3 ,]R) with compact support, for which ^| X3= o = and V^| X3= o = 
(0,0, —A 3 )\ X3=0 . A straight forward computation shows that 

Using this we see that u is an outgoing solution to the problem (12.221) if and 
only if u = e~ l ^u is an outgoing solution to the problem, 

{L A+ w, q ~ k 2 )u = e-**f in R 3 ., 

We can thus, without loss of generality, assume that A 3 = along 9M 3 , 
when showing that the solution of problem of (12.221) exists and is unique. 
We have the following result. 

Theorem 2.16. Let A and q satisfy IEE2I]) Then for any f G L 2 comp (W), 
there exists a unique outgoing solution u G Hf oc (M?_) to the problem (I2.22p . 

Proof. We shall first prove the existence. In doing so, we shall reduce the 
problem (I2.22p to all of M 3 by making use of appropriate even and odd 
extensions of the coefficients of the operator L^ q . The discussion preceding 
the Lemma 12.151 shows that we may without loss of generality, assume that 
A 3 = along <9M 3 . 

We extend the potentials A = (A±, A 2 , A 3 ) and q, and the source term / 
to the whole of M 3 . Let x := (xi,x 2 , —x 3 ). For Ai, A 2 , and q, we do even 
extensions in x 3 , i.e., 



q(x) 



Aj(x), x 3 > 0, 

q(x), x 3 < 0, 
q(x), x 3 > 0. 
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For A 3 and /, we do odd extensions in x 3 

AM = 



fa) 



A 3 (x), x 3 < 0, 
-A 3 (x), x 3 > 0, 

f(x), x 3 < 0, 
-f(x), x 3 > 0. 



Since A 3 = when x 3 = 0, we see that A G W^ p (M. 3 ,M). Furthermore, 
and/GLL^M 3 ). 

By Theorem 12.101 and Theorem 12.141 the problem 

(L A .-k 2 )u = f in IR 3 

has a unique outgoing solution u G i/ ; 2 oc (IR 3 ). 

Next we want to show that u is odd in x 3 . To that end it is convenient 
to write, 

L A - = -A - 2iA ■ V + p, p= -iV • i + i 2 + g. (2.27) 

Here one sees easily that the operators A, ^3^3 and p all preserve the parity 
in x 3 of a function that they operate on. Hence, the operator L A - — k 2 
preserves the parity in x 3 . 

Decompose u into an even and odd part with respect to x 3 , i.e. 

u = u e + u a , 

where 

Ue{x) = 7}(u(x) + U(x)), U (x) = ^(u(x) - u(x)) . 



Then 

/ = ( L A,g ~ k2 )u e + {L A - - k 2 )u , 
and using that / is odd with respect to x 3 , we conclude that 

{L A - - k 2 )u e = in R 3 . 

Now a direct computation shows that the function x 1— >■ u{x) is outgoing on 
M 3 , since w has this property. Thus, u e is outgoing, and by Theorem 12.101 
u e = 0. Hence, u is odd in x 3 . 

The Sobolev embedding theorem shows that u is continuous in M 3 , since 
u G //^(IR 3 ). Hence, u\g^3 = 0, so that u| R 3 is a solution to the half space 
Dirichlet problem fl2T22|) . 
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In order to prove uniqueness, we assume that u G H? oc (M?) is an outgoing 
solution to the problem f 1 2 . 2 2 j) with / = 0. We need to show that u = 0. To 
that end, let us consider the odd extension of u with respect to x 3 , i.e. 

m = \ U(X )\ 13 < (2^28) 
y-u(x), x 3 >0. 

Notice that since u = along 23 = 0, the function u is continuous across 
£3 = 0. 

Let us now show that u satisfies the equation, 

(L A ~ - k 2 )u = in M 3 , (2.29) 

with A and q as in the first part of the proof. Indeed, computing the first 
order partial derivatives of u, given by (I2.28|) . in the sense of distributions 
on M 3 , we obtain that 

(dju)(x), x 3 <0, 



dju(x) = J /y 3 = 1,2, 

I -(diii)(x), x 3 > 0, 

(2.30) 



-(dju)(x), x 3 >0 

(d 3 u)(x), x 3 <0, 
(d 3 u)(x), x 3 >0. 



d 3 u(x) 



Hence, we see that u G Hl oc 

One has to be more careful when computing the second order partial 
derivatives of u. For this reason, we shall give the details of the computa- 
tion below. Let ip G C^°(1R 3 ). Then denoting by (•, •) the duality between 
distributions and test functions, x' = (xi,x 2 ), and d 2 x , — d\ + d 2 , we have 

(-Am, ip) — — / uA(pdx = — / / u(x', x 3 )(dP, + d 3 )ipdx'dx 3 
Jr3 J^oo J R 2 

+ / / u(x', -x 3 ) (c^, + dl)(pdx'dx 3 
Jo Jm 2 

f P + OO P 

(d 2 ,,u)(x',x 3 )ipdx'dx 3 + / / (d 2 x , u){x' , — x 3 )(pdx'dx 3 



+ I I(x')dx', 
Jm? 
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where 



/O r+oo 
u(x', x 3 ) d\ (pdx 3 + / u(x', — x 3 ) <9 3 tpdx 3 
■00 Jo 

d 3 u(x', x 3 ) d 3 tpdx 3 - u(x', x 3 ) d 3 (p\°_ ^ 
f+00 

+ / (<9 3 u)(:r, -x 3 ) <9 3 ^dx 3 + -x 3 ) <9 3 <^|o °° 
-0 

d\ u(x', x 3 )(pdx 3 + <9 3 u(x', x 3 )(p\°_ ^ 

-00 
-00 

{d%u){x', -x 3 )ipdx 3 + (d 3 u)(z', -x 3 V|+°° 

r+oo 

dlu(x',x 3 )ipdx 3 + / (<9 3 w)(x', -x 3 )ipdx 3 . 



Hence, we have 



- (Au)(x) = r {Au)ix) > X3< °' (2.31) 
1 ,K i \(Au)(x), * 3 >0, y 1 

in the sense of distributions. Using (12.271) . (12.301) and (I2.3ip . we obtain that 

-Am(i) — 2iA(x) ■ V«(x) + {p{x) — k 2 )u(x), x 3 < 0, 



{Hi - k 



u 



(Au)(x) + 2iA(x) ■ Vu(x) + (-p(x) + k 2 )u(x), x 3 > 



We have therefore verified that u solves (I2.29p . 

Taking into account the fact that u is outgoing and applying Theorem 
|2~TU1 we finally get u = in R 3 _. 

□ 

It is now straightforward to establish the main result of this section, 
Theorem 11.11 

Proof of Theorem 11.11 Uniqueness follows from Theorem 12.161 since 
it implies that V\ — v% = for any two solutions V\ and t>2 of problem (II. 3p 
and (PP . 

It remains to check the existence of a solution. To that end, given / e 
# c 3 omp(<9M 3 _), let F E H 2 omp (R 3 ) be such that F\ dR3 = f. Theorem EH gives 
the existence of an outgoing solution v to 

(L A , q - k 2 )v = -(La,, - k 2 )F, 
v\ m s = 0. 
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Then u = v + F solves the problem (11.31) and satisfies (11.41) . The proof of 
Theorem 11.11 is complete. 



□ 
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3 The inverse problem: Proof of Theorem [L2 



The main task of this section is to prove Theorem 11.21 It will be convenient 
to set 

B_:=R 3 _nB, I := dR 3 _f] B, 

where B C R 3 is the open ball of Theorem 11.21 containing the supports of A,- 
and qj, j = 1,2. Recall that we assume that 

(OR 3 . \B) HTj ^ 0, j = 1,2. 

We can thus choose Fj, such that 

f j C Fj, F j CC dR'i \B, j = 1, 2. 

Then it follows from (11 .7ft that 

AA 1 ^(/)|f 1 = A^ iB (/)| flJ (3.1) 

for any / £ H 3 ^ 2 (dR 3 _), supp(/) C F 2 . In order to prove Theorem 11.21 we 
shall only use the data (13.11) . which turns out to be enough to determine the 
magnetic field and the electric potential. 

The gauge invariance of the DN-map plays an important role in the sequel. 
We state therefore the following results. 

Lemma 3.1. Let A £ W^I^M 3 ) and q £ L°°(W_). Then 

(i) For all if; £ C^R^R) we have 

e-^A A , q e^ = A A+ ^, q . 

(ii) There exists ip £ C 1 ' 1 (M 3 ,M) with ip\{ xs =o} = 0, for which 

and (A + W)|{« 3 =o} = (A u A 2 ,0). 
Proof, (i). A straight forward computation shows that 

e-^L A>q e^ = L A+V ^ q . (3.2) 

So that a function u solves L A+ ^^ >q u = if and only if v = e l ^u solves 
L A , q v = 0. Moreover we have 

(e-^A A>q e^)f = e-^(d n +m ■ A)(e^u) 
= id n ipu + d n u + in ■ Au 
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(ii). By Lemma \2. 151 there exists a ip E C 1 ' 1 (M. 3 _,M) with ^|{ X3=0 } = and 
V"0|{x 3 =o} = (0,0, — A3). By part (i) we have 

Aa+W = e-^A A , q (e^)f 

□ 

Remark 3.2. Notice that part (ii) of the Proposition says in other words 
that we can change a potential A to A = (A±, A2, A3 — d^ip) = (A±, A 2 ,0), 
while still retaining that Aa, q = A^ . It follows that we can always assume 
that n ■ A \ { X3 =o} — 0, where n is unit normal to the plane {x-j = 0}, without 
altering the DN-map. 

Notice also that this is the reason why we can take the DN-map as having 
the value (d n +in ■ A)u\{ X:j= o}, instead of just d n u\{ X:j= o}, even though we do 
not "know" A on the boundary. 

3.1 An integral identity 

One central step in the ideas that are used in proving uniqueness results 
for Calderon's problem, is to derive an integral equation that expresses L 2 
orthogonality between the product of two solutions ui and u 2 , and the dif- 
ference of two potentials qi and q 2 , see [35]. One shows that 

J (<?i - q2)uiu 2 = 0, 

provided that the DN-maps for q 1 and q 2 are equal. 

A similar thing will be done in this subsection, for the magnetic case. 
The integral equation, is however more involved in the case of a magnetic 
potential and cannot by itself be interpreted as an orthogonality relation. 
We will be considering the integral equation in conjunction with solutions 
that depend on a small positive parameter h. In the later sections we will 
see that in the limit h — > 0, we obtain a criterion for the curl being zero. 

We now begin deriving the integral identity. We assume that Aj,qj and 
Tj are as in Theorem 11.21 and that 

AAi,<j-i(/)|ri = A j42 , (?2 (/)|ri, 

for any / 6 Hmm P {d M. 3 _) , supp(/) C T 2 , so that (13. ip also applies. 
Let u\ G H? oc (M. 3 _) be the radiating solution to 

(L Aim - k 2 ) Ul = 0, in M 3 _, 

u l\dR 3 _ = /) 
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with / G H 3 / 2 {dRl), supp(/) C f 2 . Let v G H 2 oc (M:i) be the radiating 
solution to 

{L A2m - k 2 )v = 0, in R 3 _, 

v\dR 3 _ = /■ 

Define w := v — U\. Then 

(La 3i <b - k 2 )w = 2i(A 2 - A) ■ V«i + z'V • (A 2 - Ai)ui 

+ {A\ - Al) Ul + ( qi - q 2 )u x . (3.3) 

It follows from ( 13. ip that 

(d n + iA x ■ n)ui| fi = (d n + iA 2 ■ n)v\ ti . (3.4) 

By Remark 13.21 we may assume that A\ ■ n = A 2 ■ n = on dM. 3 _, so that 
d n w — on Ti. We also conclude from ( 13 .3p that w satisfies the equation 

(-A-k 2 )w = in R 3 _ \ IT. 

As iu|p = <9 n w|p i = 0, by unique continuation, we get that w = in M?_ \B_. 
See Theorem 14.61 and Corollary 14.71 in the appendix. Since w G H 2 oc (M. 3 _), we 
have 

w = d n w = on 9B.nl 3 .. 

Let u 2 G H 2 (B_) be a solution to (L^^ — k 2 )u 2 = in £>_. Then by Green's 
formula, we get 

(( L A 2 , q2 - k 2 )w,u 2 ) L 2 {BJ = (w, [L A2M - k 2 )u 2 ) L 2 {BJ 

- ((d n + %A 2 ■ n)w, u-^l^qbj 

+ (w, (d n + %A 2 ■ n)u 2 ) L 2 {dB _) 

= -(d n w,u 2 ) L 2^. 

Assuming that 

u 2 = on I, 

we conclude that 

(( L A 2 , q2 ~ /c 2 )w,m 2 )l2(bj = 0. 
Using equation (13. 3p we may write this as follows, 

{2i(A 2 - A) ■ (Vu)u2~ + iV ■ {A 2 - A)uiui) dx 

+ / (A\ - A\ + qi - q^uiWidx = 0. 
Jb 
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Using the fact that (A 2 — Ax) ■ n = on <9I?_ and an integration by parts, 
we get 



i / V • (A 2 — A{)uiU2 dx = — i / (A 2 — A{) ■ (Vmi«2 + u{Vu 2 )dx. 

J B_ J B_ 

Thus, we obtain that 

/ i{A 2 — A\) ■ (VU1U2 — U1VH2) dx 
Jb_ 

+ / (Al — Al + gi — q 2 )uxU2~dx — 0, 



(3.5) 



B 



where U\ G W^R 3 .) and u 2 G W 2 (B_). Here 

W^R 3 ) := {u G 4(F) I (L Al)9l - k 2 )u = in R 3 , 
supp(wi| 9ffi 3 ) C T 2 ,u radiating}, 

and 

W 2 {B_) := {u G # 2 (5_) I (L A2 ^ - k 2 )u = in «|, = 0}. 

We shall next extend the integral identity (13.51) to a richer class of solu- 
tions to the magnetic Schrodinger operators. To that end, let us introduce 
the following space of solutions, 

W X {B_) := {u G H 2 (B_) I (L Auqi - k 2 )u = in = 0}. 

The following Runge type approximation result is similar to those found in 
and 



Lemma 3.3. The space V\ := Wi(R 3 _)|b_ is dense in W\(B_) in the L 2 (B_)- 
topology. 

Proof. Suppose that V\ is not dense in Wi(B_). First notice that span(Vi) = 
V\ so that V\ is a linear subspace of L 2 (B_). Since V± is not dense in Wi(B_), 
we have a vector u G W\{B_) such that u V\. We can decompose uq 

as u = a + b, where a G V±, b G V\ and 6 / 0. Let T be the linear 
functional on L 2 (B_), defined by T(x) := proj^-±(x)/||6|| i 2, where projp^± is 

the orthogonal projection to Vi~ . Now clearly ||T(u )IU 2 = 1 an d T\ Vl = 0. 
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By the Riesz representation theorem, there is G L 2 (B_) that cor- 
responds to T. Extend g? by zero to the complement of B_ in M?_. Let 
U G H 2 oc (M. 3 _) be the incoming solution to 

(L Mm - k 2 )U = g T in R 3 _, 
U\ am s = 0. 



The existence of such a solution follows from Theorem since we can find 
a U outgoing with (L_A 1 ,q 1 — k 2 )U = g~r and U\ dR s =0. It then suffices to 

take U = U. 

Now let u G Wi(R?.). Then because T\ Vl =0 and supply) C we get 
by the Green's formula of Lemma 14.21 that 



= (u,gr)u»(a*_) = 0, (-EMi.gr ~ k )U) 

= ((L Al , qi ~ k 2 )u, U) L 2 {r3 j 
- (u, (d n + %A X ■ n)U) L 2 {m a_) 
+ ((d n + %A X ■ n)u, U) L 2 idR 3_ } 
= -{u,d n U) L2( p 2) . 

Since the boundary condition u\f can be chosen arbitrarily from C^°(F2), 
we get that d n U\f = 0. Since U\f = 0, we apply the unique continuation 

principle to conclude that U\ R a\jf — 0. As U G H 2 oc {R 3 _), we have 

U\dB_nR 3 _ — d n U\ dB _ nR 3_ = 0. 

Now applying Green's formula and doing the same calculation as above 
for Mo and £>_ instead of u yields 

(uo,9t)l 2 {bj = {uo, (Emi,3i - k 2 )U) L 2( B j 

= ((Lauq! - k 2 )u , U) L 2 {B j 
- (u , (d n + %A\ ■ n)U) L 2 {dB _) 
+ ((<9 n + iAt ■ n)u , U) L 2 {9B j 

= -(u ,d n U) L 2 {l) = 0. 

Here we have used that uq\i = 0. It follows that T(uq) = 0. This contradic- 
tion completes the proof. □ 
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Since (A 2 — A{) ■ n = on dB_, we can rewrite (13.51) in the following 
form, 

-/ miV ■ ((A 2 - At)^) dx - / i(A 2 - A x ) ■ {u x Vui) dx 
Jb Jb 



+ / (A\ — A\ + <?i - q2)u 1 u 2 dx = 0. 
Jb 



Hence, an application of Lemma 13.31 implies that the integral identity (13 . 5[) 
is valid for any U\ G W\{BS) and any u 2 G W 2 (B_). 

We summarize the discussion in this subsection in the following result. 

Proposition 3.4. Assume that Aj,qj and Tj, j = 1,2 are as in Theorem 
and that the DN-maps satisfy 

A^mir^A^/)!^, (3.6) 

for any f E H?/,l p (dm. 3 _), supp(/) C T 2 . Then 

i(A 2 — Ai) ■ (VW1M2 — U1VU2) dx 



L 



(3.7) 



+ / (A 1 - A 2 + q x - q 2 )u x u 2 dx = 0, 
Jb_ 

for Ui £ Wi(B_) and any u 2 G W 2 (B_). 



□ 



Remark 3.5. Notice that the proof of Proposition 3.4 only uses the assump- 
tion (13. ip , which follows from (13.61) . Proposition \3.4\ holds therefore also 
under the weaker assumption. 

The next step in proving Theorem 11.21 is to use the integral identity (13. 7p 
with Uj, j = 1,2, taken to be special solutions, which are called complex 
geometric optics solutions. 



3.2 Complex geometric optics solutions 

Let Oci 3 bea bounded domain with C^-boundary, and let A G W 1,0 °(n, M 3 ), 
q G L°°(f2,C). The task of this subsection is to review the construction of 
complex geometric optics solutions for the magnetic Schrodinger equation, 

L A)q u = in Q. (3.8) 
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A complex geometric optics solution to (I3.8p is a solution of the form 

u(x, C; h) = e x </ h (a(x, C; h) + r(x, C; h)), (3.9) 

where ( G C 3 , £■£ = 0, a is a smooth amplitude, r is a remainder, and h > 
is a small parameter. 

In the case when A G C 2 (f2) and g G L°°(f2), such solutions were con- 
structed in j9] using the method of Carleman estimates, and the construction 
was extended to the case of less regular potentials in [17], see also [T9] . 

Let tp(x) = a ■ x, a G 1R 3 , |ot| = 1. The fundamental role in the con- 
struction of complex geometric optics solutions is played by the following 
Carleman estimate, 

h\\u\\ HUn) < C\\e^ h h 2 L Atq e-^ h u\\ L 2 {n) , (3.10) 

valid for all u G C^°(Q) and < h < ho, which was proved in [H] and [T7] . 
Here (n) = || u IU 2 (n) + ||^Vw||x2(q). For the convenience of the reader, 

we shall present a derivation of f)3.10p in the appendix. 

Based on the estimate (I3.10p . the following solvability result was estab- 
lished in [T7J Proposition 4.3]. See also the discussion in the appendix. 

Proposition 3.6. Let A G W^°°(n, M 3 ), q G L°°(fi,C) ; a G M 3 , \a\ = 1 and 
<p(x) = a ■ x. Then there is C > and ho > such that for all h G (0, /io]> 
and any / G L 2 (Q), the equation 

e v/h h 2 L Aq e~ v/h u = f in Q, 

has a solution u G H l (Vt) with 

C 

IMIi^n) < -^II/IIl 2 (^)- 

Proof. See appendix. □ 

Our basic strategy in constructing solutions of the form (I3.9P is to write 
(JSU), as 

L c r = -L c a, (3.11) 

where := e~ x '^/ h h 2 LA :g e x '^ h - Then we first search for a suitable a, after 
which we will get r by Proposition 13.61 We must however take some care in 
choosing a and the way it depends on h, since we need later that \\r\\ H i ^ — > 
0, sufficiently fast as h — > 0. We need a also to be smooth enough. This will 
be handled as in [T7] . 
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We extend A E W 1,00 (Q, M 3 ) to a Lipschitz vector field, compactly sup- 
ported in Q, where Cl C M 3 is an open bounded set such that Q CC O. We 
consider the mollification := A * ip e E C^°(f2, M 3 ). Here e > is small 
and ip e {x) = e~ 3 i/j(x/e) is the usual mollifier with ip E C , ^°(IR 3 ), < ip < 1, 
and J ipdx = 1. We write A^ = A — A$. Notice that we have the following 
estimates for A b , 

\\A b \\ L ^ {n) = 0(e), (3.12) 
II d a A%o 0{n) = 0(e-l Q l) for all a, 

as e — t- 0. 

We shall work with a complex C = Co + Ci depending slightly on h, for 
which 

( . ( = 0, Co := a + i(3, a, (3 E S 2 , a ■ (3 = 0, (3.13) 
Co independent of h and (i = O(h), as h — >■ 0. 

By expanding the conjugated operator we write the right hand side of ( 13. lip 

as 

L ( a = (-h 2 A - 2i(-i( + hA) ■ hV - 2d • hV + h 2 A 2 

- 2ih( ■ (A a + A b ) - 2z/iCi • A - ih 2 (y ■ A) + h 2 q)a. (3.14) 

Now we want a to be such that this expression decays more rapidly than 
0(h), as h -> 0. 

Consider the operator in (I3.14p . ignoring for the time being a and its 
possible dependence on h. We would like to eliminate from this operator the 
terms that are first order in h. Notice first that Ci E O(h) and that we can 
control ||^4 b ||L oo (n) with h, if we choose e to be dependent on h. Then in an 
attempt to eliminate first order terms in h, it is natural to search for an a 
for which 

Co • Va = -i( ■ A*a, in fi. (3.15) 

We will look for a solution of the form a = e*. The above equation 
becomes then 

( . V$ = -zCo ■ A\ in ft. (3.16) 

Pick a 7 G S 2 , such that 7_L{a, /?}. 

Next we consider the above equation in coordinates y, associated with the 
basis {a, (3, 7}. Let T be the coordinate transform y = Tx := (x-a, x-(3, x-j). 
Using the chain rule and the fact that T _1 = T*, one gets thatd 

V($ o T- l )(Tx) = T[V$(x)]*. 

2 Hcrc T* is the transpose of T. 
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We therefore have that 

(1, i, 0) ■ V($ o T~ x ) (Tx) = (1, z, 0) ■ T[V$(x)]* 

= (a • V + i/3 • V)$(z) 
= Co-V$(a;). 

Equation (13.161) gives hence the (9-equation 

2 «9 2 - •($ o T- 1 )^) = -<o ■ (A» o T~ 1 )(y), (3.17) 

where <9 2 = (c^ +id y2 )/2. We will solve this using the Cauchy operator 

iV-VCs) := - [ —^—f{x - s 2 , 0))d 5l ds 2 , 

7T Si + 2S 2 



which is an inverse for the <9-operator, N := (d Vl + id V2 )/2 (see e.g. [T3] 
Theorem 1.2.2). We will need the following straight forward continuity result 
for the Cauchy operator. 

Lemma 3.7. Let r > and / G W fc '°°(]R 3 ), k > and assume that 
supp(/) C 5(0, r). Then 

1 1 AT 1 f\\w k >°°(M 3 ) — Ck\\f\\w k >°°(M. 3 ) 

for some constant C k > 0. And if / G C (M 3 ), then N~ x f G C(M 3 ). 

Proo/. See e.g. [32]. □ 



Returning to fl3~T7j) we get that $ = |A" _1 (-< • (A'or^oT. Or more 
explicitly that 

$(^,Co;e) = / cisiciss, (3.18) 

27T J R 2 Sl + 2S 2 

here T _1 (si, s 2 , 0) = S\a + s 2 f3. We have thus found a solution a = e* to 
equation (13.151) . We will choose e so that it depends on h, which implies that 
a will depend on h. In order to determine how the norm of r will depend on 
h and also for later estimates, we will need to see how || d a a||xoo depends on 
h. 

Lemma 3.8. Equation (13.151) has a solution a G C°°(Q) satisfying the esti- 
mates 

\\d a a\\ L o a(n) < C a e~ H for all a. (3.19) 
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Proof. Existence of a solution is a consequence of the considerations above. 
Therefore we need only to prove the norm estimate. 

For a = 0, Lemma 13.71 gives that ||$||.L°°(n) < C. From this it follows 
that ||e*||ioo(n) < C. For \a\ > 1 argue similarly using the estimates (I3.12p . 

□ 

We can now write the L°°(Q) norm of (I3.14p as 

||£c a IU°°(«) = II ~ h 2 L A<q a + 2ih( ■ A^a + 20 • hVa + 2ihd ■ Aa\\ L oo^y 

Using (13TT2]) . ( Eg) and the fact that Ci = 0(h) we have that 

||£c a IU°°(n) = 0(h 2 e~ 2 + he). 

Choosing e = h 1 ^ 3 , gives finally ||L^a||x,oo(n) = 0(h 4 ^ 3 ), as h — > 0. 
Finally to solve ( 13. lip for r, we rewrite it as 

e -x-Re(/h h 2 LA ^x-Re(/h( e ix-lm(/h r ^ = ^ix-JmC/h^ (3 20 ) 

If we replace e ja>Im CM r by f, then the solvability result 13. 6[ shows that we can 
find a solution f, so that a solution r to (I3.20p is given by r = e ~ iX - lm C/ h f m 

To get a norm estimate for r, notice that for the left hand side of f l3.20p 
we have 

||e^ Im ^L c a|| LO o (n) = 0(/ i 4 / 3 ), 
as h — > 0. The solvability result 13.61 gives then that 

\\r\\ HUn) = 0(h^% 

as h — > 0, which implies that ||r||#-i (q) = 0(h l l 3 ), as h — > 0. 

Thus we have obtained the following existence result for complex geomet- 
ric optics solutions. 

Proposition 3.9. Let A G W 1 ' 00 ^, M 3 ) and q G L°°(J7, C). Then for h > 
small enough, there exist solutions u G , of the equation 

La jQ u = m f2, 

t/iai are of the form 

u(xX;h) = e x </ h (a(xX;h)+r(xX;h)), 

where ( G C 3 , is of the form given by f)3.13p . a G C°°(Q) solves the equation 
f)3.15p . and where a and r satisfy the estimates 

\\d a a\\ L ~ m <C a h-\ a \/ 3 and \\r\\ HLm = 0{h 1 ' 3 ). 
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□ 

Remark 3.10. In the sequel, we need complex geometric optics solutions be- 
longing to H 2 (Q). To obtain such solutions, let Q' DD Q be a bounded domain 
with smooth boundary, and let us extend A e W 1,oc (Q, M. 3 ) and q G L°°(Q) to 
W 1,00 (Q', IR 3 ) and L°°(f2 / )- functions, respectively. By elliptic regularity, the 
complex geometric optics solutions, constructed on Q' , according to Proposi- 
tion EH belong to H 2 (tt). 

Remark 3.11. Recall that $ = ±7V _1 (-i(a + i/3) ■ (A B o T" 1 )) o T. Lemma 
3.1 implies that N~ l : Cq(Q) — > C(Q) is continuous. The estimates (I3.12p 



show that A$ — > A uniformly on Q. It follows that there is an $° ; s.t. 

||$(x,Co;/i 1/3 )-^°IU-(n)^0, h^O, 
where $° = ^N~ l (— i(a + if3) ■ (A o T -1 )) o T solves the equation 

( . V$° = -< • A in ft, (3.21) 

as h — > 0. 

Remark 3.12. We shall later use a slightly more general form for the am- 
plitude a in the CGO solutions. Namely we suppose that a = ge® , where 
g e C°°(n) ; with 

Co • V<? = 0. (3.22) 

This means that g is holomorphic in a plane spanned by a and (3. Notice 
also that by picking a = ge® , we get by (13. 15ft that 

Co ■ gV$> = -<o ■ gA\ 

in place of (I3.16p . But the $ solving (13.161) also solves the above. Hence we 
can use the same argument to obtain the $ for the above equation, as earlier. 
We thus obtain CGO solutions of the form 



u = e x</h {ge" + r 



where $ is the same as when a is of the earlier of form with no g. 

Notice also that setting a = ge* does not affect the norm estimates on a 
in Proposition ^. 9\ since g does not depend on h. 
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3.3 Recovering the magnetic field 



The aim of this section is to prove the first part of Theorem 11.21 by showing 
that the curl of the magnetic potential is determined by the DN-map. We 
use again similar notations as in Subsection 13. 1[ i.e. 



where B is an open ball in R 3 , containing the supports of the potentials Aj 
and qj, j = 1, 2. The first step in the argument will be to construct complex 
geometric optics solutions u\ and u%, belonging to the spaces Wi(B_) and 
W£ (-B-) (defined in Section |3~T1) and then to examine the limit of (13. 7p as 



For Mi E Wi(B_) and u 2 E W%{B_), we have that Uj\i = 0, j = 1,2. 
To obtain solutions that satisfy this condition, we will first choose solutions 
defined on the bigger set B = B + U / U B_. 

The parameters ( for the complex geometric optics solutions will be picked 
as follows. We will assume that 

£ ; 7i,72 E R 3 , |7i| = I72I = 1 and that {71, 72,0 is orthogonal. (3.23) 

Similarly to [31], we set 



so that Q ■ Q = 0, j = 1, 2, and (£1 + C,2)/h = i£. Here h > is a small 
semiclassical parameter. 

We need to extend the potentials Aj and qj, j = 1,2, to B + . For the 
component functions A/,i, Aj^, and qj, we do an even extension, and for 
Aj^, we do an odd extension, i.e., for j = 1, 2 we set, 



B_ := R 3 _ n B, B 



+ 



R^DB, l:=dM. 3 _nB, 




(3.24) 




x 3 < 0, 
x 3 > 0, 



£3 < 0, 
x 3 > 0, 




x 3 < 0, 
x 3 > 0, 
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where x := (xi, x 2 , — X3). By Remark 13.21 we can take Aj s\ X3= o = 0, from 
which it follows that Aj G W X% °°{B) and % G j = 1,2. 

We can now by Proposition 13.91 and Remark 13. 101 pick complex geometric 
optics solutions U\ in H 2 (B), 

of the equation (L^ - — k 2 )ui = in B, where $1 G C°°(B). By Remark 
I3~TT1 $ 1 — t- $° in the L°°-norm as h — > 0, where solves the equation 

(7i + *72)-V$? = -z(7i + i7 2 )-ii in B. (3.25) 

To obtain a function that is zero on the plane £3 = 0, we set 

U\(x) :— Ui(x) — Ui(x), xeB-UI. (3.26) 

Then it is easy to check that U\ G W\{B_). 

We can similarly pick by Proposition 13.91 and Remark 13. 10} complex ge- 
ometric optics solutions u 2 in H 2 (B), 

of the equation {L^-~- — k 2 )u 2 = in B, where $ 2 G C°°(B). By Remark 
13. llj ^2 — ^ ^2 m the -^°°-iiorm as /i — > 0, where $° solves the equation 

(-71 + ^72) • V$£ = -i(-7i + 272) • ii in B. (3.27) 

To obtain a function that is zero on the plane £3 = 0, we set 

u 2 (x) := u 2 (x) - u 2 (x), xeB^Ul. (3.28) 

Then it is easy to check that iti G W 2 *(B_). 

The next step is to substitute the complex geometric optics solutions u\ 
and u 2 , given by (I3.26|) and (I3.28p . respectively, into the integral identity 
( 13. 7p . This will be done in the Lemma bellow. We will use the abbreviations 
Pi(x) := e* 1 ^ + 77 (x) and P 2 (x) : = e* 2 ^ + r 2 (x), so that 

Ul {x) = e x<l/h P 1 (x) - e £<l//l Pi(x), 
u 2 (x) = e x<2/h P 2 (x) - e x<2/h P 2 (x). 

For future references, it will be convenient to compute the product of the 
phases that occur in the terms Uiu 2 , Vu\U2 and Mi Vw-J 

e x-Ci/h e x-( 2 /h _ e ix-£, e x-Cl/h e &-C2/h _ ^x-i 

e x<l/h e x-&/h _ e ix^ e i(0,0,-2x 3 )-( 1 /h _ e ix-t--2~fx, 3 x 3 /h (3.29) 
e -£<i/h e x-C,2/h _ e ix-£ e i(0,0,2x 3 )-d/h _ e ix-£ + +2-yi, 3 x 3 /h 

11 



where 7, = (7,- 1, 7.7,2, 7.7,3), 3 = 1,2 and 



t 2 / /z 2 |e| 2 

We restrict the choices of 71, by assuming that 

71,3 = and 72,37^0. (3.30) 

We need these conditions for the proof of the next Lemma. The first condition 
makes the above phases purely imaginary, which avoids exponential growth 
of the terms, as h . The second condition implies that |£±| — > 00 as 
h — > 0. This will be needed since, we will use the Riemann-Lebesgue Lemma 
to eliminate unwanted imaginary exponentials. 

Finally it will also be convenient to explicitly state the following norm 
estimates, which follow from Proposition 13.91 

11^11^ = 0(1), \\Ve*>\\ L ~ = 0(h^% 
llrJ^C^ 1 / 3 ), ||VrJ L2 = 0(/r 2 / 3 ), j = 1,2, 

as h — > 0. 



Lemma 3.13. If the assumptions of Proposition \3J^ hold, then 



(71 + H2) ■ I (A 2 - A 1 )e ix *e^ + *°dx = 0, (3.32) 
Jb 

where 71,72 and £ satisfy (I3.23|) and (I3.30p . 

Proof. We will prove the statement by multiplying the integral equation of 
Proposition 13.41 by h, when Ui and m 2 are given by (13.261) and f)3.28p . and 
then take the limit as h — > 0. 

We first show that for the second term in (13 .7p we have 

h [ {A\-A 2 2 + qi - q 2 )u l u 2 -> 0, (3.33) 

JB_ 

as h — > 0. Using the phase computations (13.291) we get that 

Ul u 2 = e^' x Pi(x)P 2 (x) - e tx<+ P 1 {x)P 2 (x) 
- e^-Pi(x)P 2 (x) + e l ^P l {x)P 2 {x). 

This is multiplied by an L°° function in (13.331) . Since we restricted the 
choice of 71 to make the exponents purely imaginary, we see easily using the 
estimates (13.311) that (I3.33j) holds. 
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Equation (13. 7j) multiplied by h, is thus reduced in the limit to 

lim J i{A 2 - A x ) ■ Vu x u 2 -h J i(A 2 - A{) ■ u^u^j = 0. (3.34) 

We will proceed by examining the first term. Using (I3.29P we write VuiH 2 

as 

V Ul u 2 ^(e ix <P 1 (x)P 2 J^ - e ix <+P 1 (x)P 2 JsJ) 



+ e ia ^VPi(x)P 2 (x) - e ix<+ VP l {x)P 2 



>p 2 


(x) 




x)) 


)P2 


(x) 



- e ix< -VP!(x)P 2 (x) + e ix <VPi(x)P 

where Q := Q ■ (0,0, —1), j = 1,2. The terms of the product that do not 
contain the factor 1/h, result in integrals similar to the one in ( 13 .33j) . And 
one sees similarly using estimates ( 13 . 3 1 [) that they are zero in the limit of 
(I3.34p . The first term inside the limit in ( I3.34[) is therefore reduced to 



lim / t(A 2 - A,) • (( 1 e ix <P 1 (x)P 2 (x) - C 1 e^-P 1 (£)P 2 l 



x 



-Cie iX < + P l (x)P 2 (x) - C 1 e ii <P 1 (x)P 2 (x)). 

Now we use the Riemann-Lebesgue Lemma to conclude that the terms with 
exponents containing £ + and £_ are zero in the limit. To see this, notice 
that by Remark 13.1 1[ we see that ||$i||L°°(_B_) < C, for some C > 0, when 

h is small enough. Estimates ( I3.31j) show that ||7*f Ux, 1 (s ) = (^(h 1 / 3 ). Hence 
||Pj|Ui(B_) < for some C > when h is small enough. Finally we have 
£± — > oo, as h — > 0, because of the restrictions (I3.30p . 
The first term in (I3.34j) is therefore 

lim / %{A 2 - Ax) ■ {( 1 e iX <P 1 (x)P^(xj + Q ie ix < P x {x)P^)) 

as h — > 0. The terms containing in the products of Pi and P 2 are, because 
of (13.31 j) . zero in the limit. The above limit is thus equal to 

lim f %{A 2 - A x ) ■ (£ ie **-€ e *i(*)+*^ + ^ e «-? e *i(s)+W) . 

J B 
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Finally we split the integral and do a change of variable in the second term 
and arrive at the expression 

lim f i{A 2 - At) ■ £ ie ^ e *i(*)+*^) j (3.35) 

for the first term of (13.341) . 

Returning to the second term in (13.341) . containing W1VW2. This is of 
the same form as the first one. By doing the above derivation by simply 
exchanging the roles of u\ and u^, we similarly see that the second term 
becomes 

lim - / i(A 2 -Ai) .£ e **« e *i(*)+*^T (3. 36 ) 

Now (1 — > (71 + 272) and ( 2 — > — (71 + 272), as h — > 0. Thus by using (13. 35ft 
with (I3.36p . we can rewrite (13 .34p as 

lim / i{A 2 -A x )> (de'^e* 1 * 1 )^ _ ^ e «* e *i(*)+*^J) 
h ~>° Jb 

= [ z(A 2 -A 1 )-( ll +t l2 )e tx <e ¥ ^ + '^ = 0. 
Jb 

□ 

The next Proposition shows that f!3.32|) holds even when the exponential 
function depending on i = 1, 2 is removed. 

Proposition 3.14. TTie equality (13.321) implies that 

(71 + ^l2) ■ [ (A 2 - A)e tx <dx = 0, (3.37) 

JB 

for 71, 72 and £ which satisfy (I3.23P and (I3.30p . 
Proo/ By fl3~25|) and fl3~27D we have that 

(71 + 112) ■ V($? + 1§) = -1(7! + ^72) • (ii - A 2 ) in 5. (3.38) 

Remark 13.121 furthermore implies that the amplitude e* 1 in the definition of 
Ui can be replaced by ge® 1 , if g G C°°(B) is a solution of 

(7i + ^72) • V# = in B. (3.39) 

Let ty(x) := + &2( x )- Then instead of (I3.32p we can write, 

(71 + 172) • f(A 2 - A 1 )ge ix <e^dx = 0. 
Jb 
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We conclude from (13.381) that 



(71 + il2) ■ (A 2 - Ai)ge 



+ H2) ■ (gVe ), 



and therefore, we get 




^( 7l + i l2 ) ■ Ve*dx = 0, 



(3.40) 



for all g satisfying (13.391) . 

We pick a 73, with I73I = 1, so that we obtain an orthonormal basis 
{71, 72, 73}- Let T be the coordinate transform into this basis, i.e. y = Tx = 
(x ■ 71, x ■ 72, x ■ 73). Set z = yi + iy 2 , so that d z = (d yi +i d V2 )/2 and 



Rewriting (I3.40p using this and a change of variable given by T we have 



for all g satisfying (j3.39j) . 

Notice that y ■ £ — 2/3^3, since ^ is in the y-coordinates of the form 
(0,0, £3). The above integral is therefore a Fourier transform w.r.t. £ 3 . Let 
g G C°°{TB) satisfy d 2 g = and be independent of 2/3. Then taking the 
inverse Fourier transform we write 



where T V3 := TB D and U y3 = {(2/1,^/2,2/3) : (2/1,1/2) £ M. 2 }. Notice that 
the boundary of T yz is piecewise smooth. Multiplying the above by 2i and 
using Stokes' theorem we get that 



(71 + ^72) • V = 2 d g . 










(3.41) 
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for all holomorphic functions g G C°°(T y3 ). 

Next we shall show that (I3.4ip implies that there exists a nowhere van- 
ishing holomorphic function F G C(T y , 3 ) such that 

F\ 9Ty3 =e*\ dTy3 . (3.42) 

To this end, we define F to be 

The function F is holomorphic away from dT V3 . As e* is Lipschitz, we know 
because of the Plemelj-Sokhotski-Privalov formula (see e.g. [IE]), that 

lim F{z)- lim F(z) = e* (2o) , z e dT V3 . (3.43) 

Now the function £ h-> (£ — z)" 1 is holomorphic on T y3 when z ^ T y3 . By 
choosing g(z) = C (C — i n (13.41|) . get therefore that F(z) = 0, when 
z $l T y3 . Hence, the second limit in (I3.43P vanishes, and therefore, F is 
holomorphic function on T V3 , such that (13.421) holds. 

Next we show that F is non- vanishing in T V3 . When doing so, let d T V3 be 
parametrized by z = 7(t), and N be the number of zeros of F in T V3 . Then 
by the argument principle, we get 



2ni J 1 F(z) 2ix% J Foi ( 2ni 7 e *o 7 ( 

To see that the last integral is zero, notice that this the winding number of 
the path e*° 7 . And that e^^^ is homotopic to the constant contour {1}, 
with the homotopy given by e s *( 7 ™, s G [0, 1]. 

Next, since F is a non- vanishing holomorphic function on T y3 and T y3 is 
simply connected, it admits a holomorphic logarithm. Hence, (I3.42p implies 
that 

{\OgF)\ dTyz = ^\dTy 3 - 

Because logF = \I> is continuous on dT V3 , we have by the Cauchy theorem, 

g^dz = / g log Fdz = 0, 



9T y . 3 JdT V3 



where g G C°°(T y3 ) is an arbitrary function such that <9 5 g = 0. Using Stokes' 
formula as in (13.411) allows us to write this as 



gd z ^dy 1 dy 2 = 0. 



T V3 
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Taking the Fourier transform with respect to 2/3, we get 

/ e iy <gd^dy = 0, 

JT{B) 

for all £ = (0,0, £3), £3 G R. Hence, returning back to the 2 variable, we 
obtain that 

(7i + H2) ■ / e ix< g{x)V^{x)dx = 0, 

where g G C°°(B) is such that (71 + 272) • Vg = in B. 
Using f)3.38p . we finally get 

(71 + 272) ■ f(A 2 - A 1 )g(x)e lx <dx = 0. (3.44) 
Jb 

Setting g = 1, we obtain (13.371) . 

□ 

By replacing the vector 72 by —72 in f)3.37p . we see that 

(71 - *7 2 ) ■ [ {A 2 - A x y x< dx = 0. (3.45) 

Hence, fl3T37j) and ( jS35D imply that 

7 • / (i 2 - AOe^^dsc = 0, (3.46) 

for all 7 G span{7i,7 2 } and all (Gi 3 such that (ET23I) and fET30D hold. 

In the proof of the next Proposition we see that (13.371) is actually a 
condition for having V x (Ax — A 2 ) = 0. This is therefore the last step in 
proving that the DN-map determines the curl of the magnetic potential. 

Proposition 3.15. Assume that Aj,qj and Tj, j = 1,2 are as in Theorem 
II. £1 and that the DN-maps satisfy 

A Altqi (f)\ ri = Aa 2 , q2 (f)\ri, 

for any f G Hcll P (dR 3 .), supp(/) C T 2 . Then 

V x A 1 = V x A 2 in B. (3.47) 
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Proof. Assume that £ G M 3 is not on the line L := (0, 0,t), i G R. Then the 
vectors 71 and 72 given by 

7i := (-6,6,0), 7i:=7i/|7i|, 

72:=£x7!, 72 := 7 2 /|7 2 |, (3.48) 

where £ x 7l stands for the vector cross product, satisfy (1 3 . 3 j) and (13.231) . 
Thus, for any vector £ G M 3 \ L, (I3.46J1 says that 

7-^(0 = 0, *;(£):= (3-49) 

for all 7 G span{ 7l , 72 }. Here x is the characteristic function of the set B. 
For any vector £ G M 3 , we have the following decomposition, 

«(0 = «£(0+«l(0> 

where Re^(£), Im^(() are multiples of £, and Revj_(£), Imwj_(£) are or- 
thogonal to £. Now we have Ret>j_(£), Imwj_(£) G span{ 7l , 72 }, and therefore, 
it follows from (jS3SD that v_l(0 = 0, for all (6l 3 \L 
Hence, v(£) = so that that 

e x u(o = 0, 

for all ( G R 3 \ I, and thus, everywhere, by the analyticity of the Fourier 
transform. Taking the inverse Fourier transform, we obtain (13.471) . □ 

3.4 Determining the electric potential 

In order to complete the proof of Theorem ll.2[ we need to show that the 
electric potential is also determined by the DN-map. Again we assume that 
Aj,qj and Tj, j = 1,2 are as in Theorem 11.21 and that the DN-maps satisfy 
(fi~Tj) . and hence ( l31|) . 

Since B is simply connected, it follows from the Helmholtz decomposition 
of Ax - A 2 and (jMZD that there exists V e C 1 ' 1 ^) with V = near dB 
such that 

Ai = A 2 + Vip in B. 

We extend ^ to a function of class C 1,1 on all of M 3 such that ^ = on 
M 3 \ B. Then 

Ax = A 2 + V^> in R 3 . 

In particular, ^ = on r\ U T 2 . It follows then from Lemma I3TT1 part (i) and 
(13.11) that for all / with supp(/) C T 2 , 

^1,91 (/)!?! = ^A 2 ,q 2 (/)!?! = ^ 2 +w,g 2 (/)lf 1 = ^1,92 (/)lfi- 
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We can now, by Remark l3.5l use this with Proposition ^. 41 That is we consider 
equation (13.71) . in the case A\ — A 2 . This gives 

(<7i - q-2)u{uidx = 0, (3.50) 



B 



for all Mi G Wx(B-) and u 2 G W^{B_). 

Choosing in (I3.50p u\ and u 2 as the complex geometric optics solutions, 
given by f!3.26j) and (13.281) . and letting h — )■ 0, we have 

/ (?i ~ q 2 )e ix< e^ x)+ ^dx = 0. (3.51) 
Jb 

By Remark 13. 121 in the definition (I3.26P of u\ can be replaced by ge* 1 if 
g G C°°(B) is a solution of 

(7i + ^72) ■ V<7 = in fi. 
Then ( I3.5ip can be replaced by 

/ (qi - q 2 )g(x)e lx <e*° ix)+ ^dx = 0. 
Jb 

Now (13.381) has the form, 

(71 + 272) • V($? + $f) = in 5, 

since we consider that Ax = A 2 . Thus, we can take g = e~(*° + *°) and obtain 
that 

/ (ft " q2)e lx< dx = 0, (3.52) 
Jb 

for all f G IR 3 such that there exist 71,72 G M 3 , satisfying f l3T23|) and fEOO"]) . 
Since for any ^ G M 3 not of the form £ = (0,0, £3), the vectors, given by 
(JS3HD, satisfy fl3^3|) and (13301) . we conclude that ( I532]l holds for all (el 3 
except those of the form £ = (0,0,^), and therefore, by analyticity of the 
Fourier transform, for all £ G M 3 . Hence, q\ = q 2 in This completes the 
proof of Theorem 11.21 
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4 Appendix 

4.1 Magnetic Green's formulas 

Let us first recall, following [9], the standard Green formula applied to the 
magnetic Schrodinger operator. 

Lemma 4.1. Suppose that Q C IR 3 is open and bounded, with piecewise C 1 
boundary. Let A G W 1,00 (fi, M 3 ) and g G T/ien we /iave, 

v) L 2 {n) - (u,L A> gV) L 2 {n) 

= (u,(d n + iA-n)v) L 2 {dQ) - {{d n + iA ■ n)u, v) L 2 {m) , 

for all u,v G H 1 {Q), with Am, Af G L 2 {Q), where n is the exterior unit 
normal to dfl. 

We shall also need a version of the above result where Q is replaced 
by ffi 3 . We shall then need to put some restrictions on v and u, because 
M 3 is unbounded. To this end we assume that u and v are solutions to 
the Helmholtz equation outside some compact set, that obey some form of 
radiation condition. To be precise, let A G W^ p {R 3 _, M 3 ), q G L^ mp (M 3 ), 
and let u G H 2 QC (M?_) be such that 

{L Aq - k 2 )u = in K 3 _, 

supp(w| 9R 3 ) is compact, and u is outgoing. Assume also that v G H 2 oc {R 3 _) 
satisfies 

(L A ^-k 2 )vEL 2 comp (W), 
supp(f | S1R 3 ) is compact, and v is incoming. 
Lemma 4.2. With u and v as above, we have 

((L A , q - k 2 )u, v) L 2 {R 3j - 0, (L A g - k 2 )v ) L 2 (R 3j 
= («, (<9„ + L4 • n)u)i3( SR 3j - {{d n + iA ■ n)u, u) L a( 8R 3j. 

Proof. Let .Br := {x G M 3 | |x| < i?} be an open ball in IR 3 of radius R, and 
choose i? > large enough so that 

supp(A),supp(g) C B R . 

Set Q = R 3 n Br. By Lemma [4. 1[ we know that 

{{L A>q - k 2 )u,v) L 2 {n) - {u, {L A g - k 2 )v) L 2 (n) 
= {u, {d n + iA ■ n)v) L 2 {m) - {{d n + iA ■ n)u, v) L 2 [m) . 
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Thus, to obtain (14.1 p we need to show that 



{ud n v - (d n u)v)dS R 0, i? ^ oo. (4.2) 

Let us rewrite the left hand side of the above as follows, 

{d n v — ikv)udS r — I {d n u — ikufvdS r. 



L 

We show that first term goes to zero as R — > oo. The second term can be 
handled in the same way. Applying Cauchy-Schwarz gives 

/ {d n v — ikv)udS r < / \d n v — ikv\ 2 dSR / \u\ 2 dSR. 

JdBnnR 3 JdBnnM 3 JdBnnR 3 



Here the first integral goes to zero, since d n v — ikv = d n v+ikv and \d n v+ikv | 2 
is o(l/r 2 ) as r = \x\ —> oo, since v is incoming. The second integral is 
bounded, because of Lemma 12.61 We conclude that (14. 2 p holds. □ 



4.2 Carleman estimates and solvability 

Let Q C M™, n > 2, be a bounded domain with C°° boundary, and let 
(p(x) = a ■ x with a G M n , |a| = 1. Consider the conjugated operator 

L 9 := e^ h h 2 L Aq e-^ h . 

Note that depends on A, q and h. In the beginning of this subsection we 
shall establish the following Carleman estimate, where we write 

1 1 U I \ 2 HlJp.) = \\ U \\l 2 (Q,) + II^^IIl 2 ^)- 

Theorem 4.3. Let A e W 1 ' 00 ^, C n ) and q E L°°(Q,C). Then there exist 
C > and h > such that for all u G C^(Q), we have 

h\\u\\H^(U) < C\\L v u\\ L 2 {n) , (4.3) 

when < h < ho. 

Proof. In what follows, the L 2 (f2)-norm is abbreviated to || ■ ||. Let e > 0. 
Define 

ip e {x) = a ■ x + -(a ■ x) 2 . 
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Denote by L 0iipe := —e ip ^ h h 2 ^e v>e ^ h . This can be decomposed as L ^ t = 
A e + iB e , where 

A e := -h 2 A - (1 + ea • x) 2 , 

B e := —2h(l + ea ■ x)a ■ iV — ieh. 

A direct calculation gives that 

II L ^u\\ 2 = \\Au\\ 2 + ||£? £ u|| 2 + i([A e: B e ]u, u), u G C °°(fi). (4.4) 

Another straightforward calculation shows that the commutator can be writ- 
ten as 

3 

i[A e , B e ] = -Aeh 3 ^ aja k djd k + Ahe(l + ea ■ x) 2 . 
j,k=i 

The first term is an operator with a positive semi-definite semiclassical sym- 
bol of the form 4eh(a ■ £) 2 > 0. The inner product in (14. 4p . with this part is 
therefore non-negative. We can hence drop it and estimate f)4.4p as follows, 

||£o,^w|| 2 > ||-4ew|| 2 + 4/ie||(l + ea ■ x)u\\ 2 

> \\A e u\\ 2 + he\\u\\ 2 , (4.5) 

for all e > such that e\a ■ x\ < 1/2, iGfl. 

The next step is to obtain a similar estimate for the first order pertur- 
bation L Ve := e^^h^LA^e-'^^ ■ We decompose this as L 9i = L 0iipt + Q e , 
where 

Q t ■= e ^ /h h 2 {-2iA ■ V - zV ■ A + A 2 + q)e~^l h . 

This can be estimated as 

||Q e u|| 2 = /i 4 || - 2iA ■ Vu + 2iA ■ V^u - iV ■ Au + A 2 u + qu\\ 2 

<C(h A \\Vu\\ 2 + h 2 \\u\\ 2 ), 

when h is sufficiently small. For L V€ , we get 

\\L , V M\ 2 < C{\\L V M\ 2 + ^ 4 ||Vu|| 2 + h 2 \\u\\ 2 ), (4.6) 

when h is sufficiently small. Next we rewrite the gradient term and use the 
Cauchy inequality to obtain 

h 2 \\Vu\\ 2 = (-h 2 Au,u) 

= (A e u, u) + || (1 + ea ■ x)w|| 2 

<^P £ n|| 2 + ^||u|| 2 + 4||n|| 2 , (4.7) 
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where 5 > 0. Combining this with (14. 5p and (14 .6p . gives 
\\A e uf + he\\u\\ 2 < c( \\L Ve u\\ 2 + h 2 \\u\\ 2 



+ 6 h 2 \\A e u\\ 2 + d + ^h 2 \\u\\ 2 
2 o 

Choose 8 — 2/C. Rearranging the above gives 

(1 - h 2 )\\A e u\\ 2 + h{e - (C 2 + 5C)h)\\u\\ 2 < C\\L v u\\ 2 - 

We may assume that h < 1/2. The first term is then larger than /i 2 ||A e M|| 2 . 
Next we pick e so that (e - (C 2 + 5(7) /i) = 6h, i.e. e = Mh, M > is fixed. 
This gives 

h 2 \\A e u\\ 2 + Qh 2 \\u\\ 2 < C\\L^u\\ 2 . 
Using (14. 7p with 5 = 2, gives then 

/i 4 ||Vwf + h 2 \\u\\ 2 < C\\L^u\\ 2 . (4.8) 
Written in another way we have 

h 2 \\u\\ 2 Hl <C\\L v ,u\\h, 

scl 

which is almost the desired estimate. 

To finish the proof we need to replace L lfie in (14.81) by L v . To this end 
let g := e(a ■ x) 2 /2 and let u = e g/h u. Notice that g/h = M(a ■ x) 2 /2 is 
independent of h. Estimate (14.81) gives now 



h 4 \\Vu\\ 2 + h 2 \\u\\ 2 < C\\e 9/h L v u\\ 2 < C'\\L 



r. II 2 



for some constants C, C > 0. To obtain (14.81) with Lp, we need only to show 
that 

h 4 \\Vu\\ 2 + h 2 \\u\\ 2 < C(/i 4 ||Vnf + h 2 \\u\\ 2 ) (4.9) 

for some constant C > 0. Using the triangle and Cauchy inequalities we see 
that 

\\e 9/h Vu\\ 2 < 2\\ V{e 9/h )u + e 9/h Vu\\ 2 + 2|| V{g /h)e 9/h u\\ 2 

< 2||V(e 9/ ^)|| 2 + C/i- 2 ||e^|| 2 , 

for some some constant C > 0. Hence, 

h^WVuf < C(h 4 \\Vu\\ 2 + h 2 \\u\\ 2 ) 1 

for some constant C > 0, which shows that (14. 9 p holds. The proof is com- 
plete. 

□ 
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Let 

imi^) = ^r 3 / (i + h 2 \er\m\ 2 dt * e r. 

We have the following consequence of Theorem 14.31 see also [9] . 

Corollary 4.4. Let A G W 1 ' 00 ^, C n ) and q G L°°(fi,C). Then there exist 
C > and ho > snc/i £/ia£ /or a// it G C^°(f2), we /lave 



folMU^n) < C\\L v u\\ H -i^ n) , (4.10) 

w/ien < h < h . 

The formal adjoint of L v is given by 

and Corollary 14.41 still holds for L* 

Next we prove the following solvability result, see also [TTJ Proposition 
4.3]. 

Proposition 4.5. Let A G W 1,0O (fi,C n ), g G L°°(fi,C) ; a G M n ; \a\ = 1 
and = a-x. Then there is C > and /io > such that for all h G (0, ho], 
and any f G L 2 {Vt), the equation 

LipU = f in Q, 

has a solution u G H 1 ^) with 

C 

\\ u \\h^ c1 (Q) < ^ll/IUw 

Proof. The Carleman estimate of Corollary 14.41 applied to L* shows that L* 
is injective on C£°(f2). This lets us define the functional T : L*C£°(fi) — > C, 
given by 

T(w) := ((L;)^wj) L2(Rny 

The set Dom(T) C C if~ 1 1 (M n ) is a linear subspace. The linear 

functional T is moreover bounded, since by Corollary 14.41 



\t(w)\ = \((l;)- 1 wJ) 



L 2 (Q)I 



< — ll/l|L 2 (n)lkl|ff-i( R n), C > 0. 
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The Hahn-Banach theorem allows us to extend T, without increasing its 
norm, to an operator T : (R n ) — > C. By the Riesz representation theorem 
there exists r e H~l(M. n ) such that 

fH = Kr)_x, iov w E H- c l(R n ), 

scl 

and 

IH| fli i (R -) = ||r||<||T||< ^11/11^(0,. (4.12) 

Here we have used (14. lip . 
Furthermore, we have 

(w,r) H -i = (w,u) L 2, 

scl 

with 

u = jT\1 + h 2 t 2 Y l Fr e Hl cl (R n ). 

Here J 7 is the Fourier transformation on R™. 

We now show that u solves L^u — f in the weak sense in Q. For every 
ip G C^°(n) we have 

= ((L^L^,f) L2{Rn) 

Hence a is a weak solution. 

To obtain the norm estimate, we observe that ( K ») 
and use f)4.12p . The proof is complete. 



□ 



4.3 The unique continuation principle 

In this work we make heavy use of the so called unique continuation principle. 
The unique continuation principle can be seen as an extension of the familiar 
property that an analytic function that is zero on some open set is identically 
zero. 

Let Q C M n be an open connected set, and let 

n 

Pu = o-ij(x) didjU + bi{x) diU + c{x)u. 

i,j=l i 
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Here G C 1 (r2) are real- valued, = ctji, and there is C > so that 

n 

J2 OiiWtej > Ciei 2 , xefi, ^1". 

Furthermore, 6j G L°°(f2, C) and c G C). We have the following result, 

see [5] and [23] . 

Theorem 4.6. Lei it G ifj^fi) &e sitc/i £/ia£ Pu = in ft. Let oj <Z Vt be 
open non-empty. If u = onw, £/ien u vanishes identically in Q. 

Corollary 4.7. Assume that dQ is of class C 2 . Let T G dQ be open non- 
empty. Let u G H 2 (Q) be such that Pu = in Q. Assume that 

u = B v u = on r. 

Here B u u is the conormal derivative of u, given by 

n 

B v u = Vi(aijdju)\ dn G H^idQ). 

»>J=1 

Then u vanishes identically in Q. 
4.4 Rellich's lemma 

Rellich's lemma is a fundamental result in the scattering theory of the Helmholtz 
equation, see see e.g. [7j. 

Proposition 4.8. Let k > and let u G V{E?) satisfy the Helmholtz equa- 
tion (—A — k 2 )u = outside a ball B in M 3 . Assume that 

lim / \u\ 2 dS^0. 

R ^°°J\x\=R 

Then u = in R 3 \ B. 
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